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Abstract. The gradients in line strength indices such as 
Mg2 and (Fe) (Carollo & Danziger 1994a, b; Carollo et al. 
1993), and their increase toward massive system (Worthey 
et al. 1992, Faber et al. 1992, Davies et al. 1993, Weiss 
et al. 1995) observed in elliptical galaxies are customar- 
ily considered as indicators of gradients in the chemical 
abundances of Mg and Fe and arguments are given for 
an enhancement of Mg (a-elements in general) with re- 
spect to Fe toward the center of these galaxies or going 
from dwarf to massive ellipticals. In this paper, we present 
a detailed analysis of the whole problem aimed at under- 
standing how the indices Mg^ and (Fe) depend on the 
the chemical abundances and age of the stellar mix of a 
galaxy. We find that Mg2 and (Fe) do not simply cor- 
relate with the abundances of Mg and Fe and the ratio 
Mg I Fe, but depend also on the underlaying distribution 
of stars in different metallicity bins N(Z), or equivalently 
on the past history of star formation. In practice, inferring 
the abundance of Mg and Fe and their ratio is hampered 
by the need of some information about N(Z). Finally, the 
observational gradients in Mg2 and (Fe) do not automat- 
ically imply gradients in chemical abundances. 

Key words: Galaxies: ellipticals - Galaxies: evolution - 
Galaxies: stellar content - Galaxies: gradients 



1. Introduction 

The line strength indices Mg^, (Fe), Hp, etc. and their 
gradients are customarily used to infer the age and metal- 
licity and their variations across galaxies. 
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Furthermore, in elliptical galaxies, the gradients in 
M<72 and (Fe) have different slopes (Fisher et al. 
1995,1996, Carollo & Danziger 1994a,b; Carollo et al. 
1993, Davies et al. 1993), which is used as a clue to ar- 
guing that Mg (a-elements in general) is enhanced with 
respect to Fe toward the center. 

The bottom line to infer from Mg^ and (Fe) an en- 
hancement in a-elements rests on the notion that these 
two indices depend on age and the abundances of Mg 
and Fe, and that age and abundance effects can some- 
how be disentangled. If this is the case, the implications 
are of paramount importance. It is worth recalling that 
according to the current nucleosynthesis scenario Fe is 
mainly produced by Type la supernovae (accreting white 
dwarfs in binary systems in the most popular scheme) 
and in smaller quantities by Type II supernovae. In con- 
trast, only Type II supernovae contribute to oxygen and 
a-elements. Furthermore as the mean lifetime of a binary 
system (Type la progenitors) is > 1 Gyr, the contami- 
nation by Type I supernovae occurs later as compared to 
Type II supernovae. Finally, we expect the iron abundance 
[Fe/H] and the [a/Fe] ratios to increase and decrease, re- 
spectively, as the galaxy ages. In standard models of galac- 
tic chemical evolution, i.e. constant initial mass function 
and supernova driven galactic winds (cf. Matteucci 1997 
for a comprehensive review of the subject), this means that 
to obtain a galaxy (or region of it) enhanced in a-elements 
the time scale of star formation over there must be shorter 
than about 1 Gyr. This is a very demanding constraint on 
models of galaxy formation and evolution. It is worth clar- 
ifying that such a conclusion is largely independent of the 
IMF and galactic wind model in usage, even if some ad 
hoc combinations of IMF and/or galactic wind model can 
be found in which enhancement of a-elements is possible 
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irrespective of the argument about the time scale of star 
formation. The reader is referred to the review article by 
Matteucci (1997) and the recent study by Chiosi ct al. 
(1997) on an unconventional IMF for more details. 

In this paper, we address the question as to what ex- 
tent the indices Mg 2 and (Fe) (and their gradients) de- 
pend on age and chemical abundances, paying particular 
attention to the complex environment of a galaxy in which 
stars of many ages and chemical compositions are present. 
In other words, we seek to clarify how the past history of 
star formation, which manifests itself in the extant rela- 
tive distribution of stars per mctallicity bin [thereinafter 
the partition function N(Z)\, affects the correspondence 
between indices, ages and abundances. 

To this aim, we will utilize the spherical models of el- 
liptical galaxies developed by Tantalo et al. (1997), which 
take into account gradients in mass density, star forma- 
tion rate, and chemical abundances (Section 2). With the 
aid of these models, we predict how the gradients in Mg 
and Fe (and their ratio) translate into gradients in Mg 2 
and (Fe), and check whether a gradient in Mg 2 steeper 
than a gradient in (Fe) implies an enhancement of Mg 
with respect to Fe toward the center of these galaxies. 
We anticipate here that, while these models are indeed 
able to match many key properties of elliptical galaxies, 
including the gradients in broad band colors (see below), 
they lead to contradictory results as far as the gradients 
in line strength indices are concerned (Section 3). To un- 
derstand the physical cause of this odd behaviour of the 
models, we check the calibration in use and the response 
of Mg 2 and (Fe) to chemistry (Section 4). The reason of 
the contradiction resides in the dependence of the indices 
in question on the existing N(Z), i.e. the past history of 
star formation. It will be shown that gradients in the in- 
dices Mg 2 and (Fe) do not automatically correspond to 
gradients in the Mg and Fe abundances, and their ratios 
in particular (Section 4). In order to cast more light on 
this topic, in a way independent of the particular model 
of galactic evolution, we derive the above indices for single 
stellar populations (SSP) of different metallicity and age, 
and look at the possible combinations of these two param- 
eters leading to the same values for the indices (Section 
5). Finally, some concluding remarks are drawn in Section 
6. 



Tq it matches the radial density distribution derived by 
Young (1976) for spherical systems. The density profile 
of the dark-matter is taken from Bertin et al. (1992) and 
Saglia et al. (1992) however adapted to the Young for- 
malism for the sake of internal consistency. The mass of 
dark matter is taken in fixed proportions with respect to 
the luminous one, Md = 9 x Ml (all models below are 
for 9 = 5), Dark matter will only affect the gravitational 
potential. 

Given the total present day luminous mass Ml of the 
galaxy (thereinafter in units of 10 12 x M Q and shortly in- 
dicated by Ml,ii) and its effective radius Rl,c, the model 
is divided into a number of spherical shells, with proper 
spacing in radius and mass. The luminous mass of each 
shell is written as 

AM LiS = p L (r) x AV(r), 

where AV(r) is the volume of the shell and p L (r) is the 
mean density of baryonic material, respectively, The in- 
ner and outer radius of each shell are chosen in such a 
way that, using the Young (1976) tabulations of the mean 
density as a function of the fractionary radius s = t/Rl.s, 
the mass AMl.s amounts to about 5% of the total mass 
Ml,i2- The radial variation of the gravitational potential 
of dark and luminous mass can be easily derived from the 
above description. The reader is referred to Tantalo et al. 
(1997) for all other details. 

Accretion rate. The collapsing rate of luminous mate- 
rial (gas) is expressed as 



dp L (r, t) 



dt 



= Pw(r) x exp{ i —) 

inf T V> 



(1) 



where r(r) is the local time scale of gas accretion to be 
discussed below, and PLo( r ) is fixed by imposing that at 
the present-day age of the galaxy Tq the density of lumi- 
nous material in each shell has grown to the value given 
by the Young profile. 



Star formation rate. This follows the standard Schmidt 



law 



dp g {r, t) 



dt 



v{r)p„{r,t) 



(2) 



sf 



2. The reference model 

Sketch of the model. Elliptical galaxies are supposed to 
be made of baryonic and dark matter both with spheri- 
cal distributions but different density profiles. While dark 
matter is assumed to have remained constant in time, the 
baryonic material (initially in form of primeval gas) is sup- 
posed to have accreted at a suitable rate onto the poten- 
tial well of the former. The rate at which the density of 
baryonic material grows with time at any galacto-centric 
distance is chosen in such a way that at the galaxy age 



where p g (r, t) is the mean local gas density and v{r) is the 
specific efficiency of star formation. 

Accretion time scale r(r). A successful description of 
the gas accretion phase is possible adapting to galaxies the 
radial velocity law describing the final collapse of the core 
in a massive star (Bethe 1990), i.e. free-fall [v(r) ccr^s] 
in all regions external to a certain value of the radius r* 
and homologous collapse inside [v(r) cx r]. This picture is 
also confirmed by Tree-SPH dynamical models of elliptical 
galaxies (cf. Carraro et al. 1997 and references therein). 
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Let us cast the problem in a general fashion by expressing 
the velocity v(r) as 

v(r) = C\ x r a for r < r* 

v(r) = C2 x r - ' 3 for r > r* 

(where ci, C2, a and /3 are suitable constants), and the 
time scale of accretion as 

r(r) cx 



as a measure of the net efficiency of star formation. 



v(r) 



In the models below we adopt a = 2 (as suggested by the 
Tree-SPH calculations) and j3 = 0.5 (as indicated by the 
core collapse analogy). The determination of the constants 
ci and C2 is not strictly required as long as we seek for 
scaling relationships. 

The time scale of gas accretion can be written as pro- 
portional to some arbitrary time scale, modulated by a 
correction term arising from the scaling law for the ra- 
dial velocity. For the time scale base-line we can take the 
free-fall time scale tff referred to the whole system, 



r(r) 



tff x — 
r 



T(r)=t ff x ( — ) 



3/2 



if r < r* 



if r > r* 



(3) 



(4) 



For the free-fall time scale tff we make use of the relation 
by Arimoto & Yoshii (1987) 



t ff = 0.0697 x Mlif 



Gyr. 



(5) 



Finally, we take r* — \Rh. e . Other choices are obviously 
possible without changing the overall results of this study. 

Specific efficiency of star formation. In order to de- 
rive the specific efficiency of star formation v{r) we utilize 
the simple scale relations developed by Arimoto & Yoshii 
(1987) however adapted to the density formalism. At the 
typical galactic densities (10~ 22 - 10 -24 g cm~ 3 ) and con- 
sidering hydrogen as the dominant coolant (Silk 1977) the 
critical Jeans length is much smaller than the galactic ra- 
dius, therefore the galaxy gas can be considered as made 
of many cloud lets whose radius is as large as the Jeans 
scale. If these clouds collapse nearly isothermal without 
suffering from mutual collisions, they will proceed through 
subsequent fragmentation processes till opaque small sub- 
units (stars) will eventually be formed. In such a case the 
stars are formed on the the free-fall time scale. In contrast, 
if mutual collisions occur, they will be supersonic giving 
origin to layers of highly cooled and compressed material, 
the Jeans scale falls below the thickness of the compressed 
layer and fragmentation occurs on the free-fall time scale 
of the high density layers, and finally the whole star form- 
ing process is driven by the collision time scale. On the 
basis of these considerations, we take the ratio 
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*// x 1 




2 3 

Age(Gyr) 



Fig. 1. Temporal evolution of four abundance ratios: [Fe/H] 
(solid line), [Mg/Fe] (dotted line), [O/Fe] (dashed line), and 
[Si/Fe] (long-dashed line). The dot-dashed line shows the ratio 
Z/Zq as a function of time. The abundance ratios are in the 
standard notation 



Let us express v(r) as the product of a suitable yet ar- 
bitrary specific efficiency v* referred to the whole galaxy 
and a dimensionless quantity F(r) describing as the ratio 
of eq. (^J) varies with the galacto-centric distance. An ob- 
vious expression for F(r) is the ratio (^|) normalized to its 
central value. 

According to Arimoto and Yoshii (1987) the mean col- 
lision time scale referred to the whole galaxy can be ex- 
pressed as 



0.0072 x Ml\ 2 



Gyr 



(7) 



With aid of this and the relation for the the free-fall 
time scale above we can first calculate v* 



tff X t r 



(8) 



J gal 



Extending by analogy the definition of free-fall and 
collision time scale to each individual region, we get 



col 



(6) F (r) = (^yx 



p„(r c ,T G ) 



P 9 0,7g) 



(9) 
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where p„(r, Tg) is the mean gas density within the region 
of mid radius r and r c is the mid radius of the innermost 
sphere. 

In principle, the exponent 7 could be derived from 
combining the mass dependence of tff and t co i, i.e. 7 ~ 
0.2. However, the many models calculated by Tantalo et 
al. (1997) show that in order to recover the observational 
gradients in broad band colors (and other properties of el- 
liptical galaxies), the efficiency of star formation (i.e. F(r) 
in our formulation) must vary with the radial distance 
more strongly than predicted by 7 = 0.2. The following 
relation is found to give good results 

7 - 0.98 x M° L f 2 (10) 

Finally, the total expression for v(r) is 

v(r) = x — x - JL - — — — Gyr (11) 

V ' UffXtcoi\ gal \rJ [p g (r,T G )\ K ' 

Remark. Before proceeding further, we like to briefly 
comment on the apparent complexity of the model 
adopted to perform the analysis. First of all, the model 
has to be considered as gross tool to understand how gra- 
dients in star formation would affect gradients in metal- 
licity and photometric properties. Secondly, the analysis 
itself is almost model-independent as what matters here is 
to make use of a reasonably grounded formulation able to 
predict gradients in chemical abundances and in narrow 
band indices and look and the mutual correlation between 
them. 

Galactic winds. The models allow for galactic winds 
triggered by the energy deposit from Type I and II su- 
pernova explosions and stellar winds from massive stars. 
The formalism in usage here is the same as in Bressan et 
al. (1994) and Tantalo et al. (1996) but for two important 
details: first only a fraction (77 = 0.3) of the kinetic energy 
by stellar winds is let thermalize, and second the cool- 
ing law for supernova remnants is the same as in Gibson 
(1994, 1996). When the total thermal energy of the gas in 
each region exceeds the local gravitational potential, gas 
is supposed to escape from the galaxy and star formation 
to stop over there. 

It worth noticing for the sake of clarity that in presence 
of galactic winds, the local density and total mass of bary- 
onic material in turn can never reach to the asymptotic 
values pl(t,Tg) and M^ia, respectively. The discussion 
by Tantalo et al. (1997) of this topic clarifies, however, 
that the final results of the models are not too severely 
affected by contradiction between the initial hypothesis 
(models constrained to match the asymptotic mass) and 
the real value of this reached in the course of evolution. 

Chemical and photometric evolution. The chemical 
evolution of elemental species is governed by the same set 
of equations as in Tantalo et al. (1996) however adapted 
to the density formalism and improved as far as the ejecta 



and the contribution from Type la and Type II supcrnovae 
are concerned according to the revision made by Portinari 
et al. (1997) and Tantalo et al. (1997) to whom we refer. 
Finally, the line strength indices Mg 2 and (Fe) have been 
calculated adopting the calibrations by Worthey (1992) 
and Worthey et al. (1994) as a function of [Fe/H], T cff 
and gravity of the stars. 

Mass-radius relationship. The final step is to adopt a 
relationship between R^,e and Ml so that assigned the to- 
tal baryonic mass, the effective radius and all other quan- 
tities in turn are known. For the purposes of this study 
and limited to the case of Hq = 50 Km sec -1 Mpc~ x , we 
derive from the data of Carollo et al. (1993) and Goud- 
frooij et al. (1994) the following relation 

R L , e = 17.13 x A/0 jf (12) 

where Rl.c is in kpc. 

Main results. This simple modelling of the distribu- 
tion of density and hence mass in a spherical system al- 
low us to describe the gradients in star formation, chem- 
ical abundances, ages, and photometric properties. These 
models indeed are able to reproduce (i) the slope of colour- 
magnitude relation by Bower et al. (1992a, b); (ii) the UV 
excess as measured by the colour (1550-V) by Burstein et 
al. (1988); (iii) the mass to blue luminosity ratio (M/Lb)© 
of elliptical galaxies. See Tantalo et al. (1997) for all other 
details. 

For the purposes of the discussion below, in Table 1 we 
present the basic data for the central core (r = 0.05i?L.e) 
and the first shell (0.05 x R L e < r < 0.15 x i?L, e ) in a 
typical galaxy with total luminous mass 3A/£,.i2- 

The content of Table 1 is as follows: AM^s is the 
asymptotic luminous mass of the region in units of 
10 12 M Q ; v is the local efficiency of the SFR; r is the local 
time scale of gas accretion in Gyr; t gw is the age in Gyr at 
which energy heating by supernova explosions and stellar 
winds exceeds the local gravitational energy; Z max and 
(Z) are the maximum and mean metallicity, respectively; 
G(t) and S{t) are the local fractionary masses of gas and 
stars gas, respectively, both normalized to the asymptotic 
mass AMl,s] ^enh is the percentage of stars showing a- 
enhancement that are present in the region (see below). 
Finally, in Fig. [j] we show the temporal evolution of the 
abundances of a few elements in the central core of the 
model. No detail for the remaining regions is given here 
but for the gradients in broad band colors and narrow 
band indices to be presented in Fig. || below. 

3. Casting the gradient contradiction 

In Fig. U we compare the theoretical and observational 
gradients in broad-band colors (left panel) and line 
strength indices (right panel) for our proto-type model. 
The data are from Carollo & Danziger (1994a,b). 
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Table 1. Basic features for the central core and first shell of 
the reference model with 3Ml,i2 



Parameter 


Core 


I s * shell 


AM L , S 


0.146 


0.150 


V 


7.1 


50.0 


T 


0.74 


0.29 


tgu 


5.12 


0.79 


7 

^max 






(Z) 


0.0365 


0.0286 


G(t) 


0.004 


0.010 


S(t) 


0.845 


0.874 


N en h 


45.7% 


53.9% 



Remarkably, while the model matches the gradient in 
broad band colors, it fails as far as the line strength indices 
Af 92 and {Fe} and their gradients are concerned. Other 
cases can be found in the Carollo & Danziger (1994a, b) 
list (they are not shown here for the sake of brevity) , in 
which either the broad band colors or the line strength 
indices are matched, but the simultaneous fit of the two 
sets of data is not possible. 

The obvious attitude toward these matters would be 
to ascribe the above failure to inadequacy of the models 
(which may certainly be the case) and thus quit the sub- 
ject. This is a point of embarrassment because there is no 
obvious explanation as to why models that successfully 
reproduce many of the observed properties of elliptical 
galaxies (cf. Bressan et al. 1994, 1996; Tantalo et al. 1996; 
1997) fail to match the line strength indices. 

In addition to this, and even more relevant to the aims 
of the present study, a point of contradiction between 
chemical structure and line strength indices is soon evi- 
dent. The problem is cast as follows. The theoretical gradi- 
ents [AlnMg 2 /AR] ~ -0.13 and [Aln(Fe)/AR] ~ -0.11 
(with R the galactocentric distance in kpc) are nearly 
identical. 

This is a surprising result because looking at the data 
of Table 1 the duration of star formation was much longer 
in the central core than in the outer shell (the same trend 
holds for all remaining shells not considered here), which 
implies that the stars in the core are on the average less 
a-enhanced than in the more external regions (cf. the 
temporal evolution of the elemental abundances shown in 
Fig. , whereas the gradients we have obtained seem to 
indicate a nearly constant ratio [Mg/Fe\. 

To single out the reason of the contradiction we look 
at the variation of the abundance ratios [Fe/H], [C/Fe], 



[O/Fe], and [Mg/Fe] (with respect to the solar value) as 
a function of the metallicity and time, and the present- 
day partition function N(Z). This allows us to evaluate 
the fraction of living stars with metallicity above any par- 
ticular value and with abundance ratios above or below 
the solar value. The relationships in question are shown 
in the two panels of Fig. || (the left panel is for the central 
core; the right panel is for the 1 st more external shell). 

In addition to this, we also look at the current age 
of the stellar population stored in every metallicity bin 
(we remind the reader that the metallicity in this model 
in a monotonic increasing function of the age, cf. Fig. |]). 
The top axis of Fig. || shows the correspondence between 
metallicity and birth-time of the stellar content in each 
metallicity bin, shortly named single stellar population 
(SSP). The SSP birth-time is t = T G - T S sp, where T G 
is the present-day galaxy age and T$s p is the current age 
of the SSP. 

From this diagram we learn that the external shell is 
truly richer in a-enhanced stars (~ 53.9% of the total) 
than the central core (~ 45.7% of the total). The percent- 
ages N en h are given in Table 1. This confirms our expec- 
tation that the models should predict gradients in line 
strength indices consistent with the gradients in abun- 
dances. 

Which is the reason for such unexpected contradic- 
tion? 

One may argue that the above disagreement results 
either from the adoption of calibrations, such as those by 
Worthey (1992) and Worthey et al. (1994), which include 
the dependence on [Fe/H], T e fj, and gravity but neglect 
the effect of enhancing the a-elements, or the particular 
model in usage. 



4. Changing calibrations and chemistry 

To answer the question posed in the previous section, 
first we adopt a different calibration in which the effect 
of [Mg/Fe] is explicitly taken into account, and second 
we discuss different, ad hoc designed, galactic models in 
which different levels of enhancement in a-elements are let 
occur by artificially changing the history of star formation. 



4--1. Calibrations containing [Mg/Fe] 

Many studies have emphasized that line strength indices 
depend not only on the stellar parameters T e // and grav- 
ity, but also on the chemical abundances (Barbuy 1994, 
Idiart et al. 1995, Weiss et al. 1995, Borges et al. 1995). 

We start pointing out that in presence of a certain 
degree of enhancement in a-elements one has to suit- 
ably modify relationship between the total metallicity Z 
and the iron content [Fe/H]. Using the pattern of abun- 
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Age(Gyr) 



Table 2. Basic properties of the central region of the test mod- 
els with 3Ml,i2- Model-A: late galactic wind and no enhance- 
ment of a-elements. Model-B: early galactic wind and enhance- 
ment of a-elements. Model-C: recent burst of star formation, 
and galactic wind, and strong enhancement of a-elements 



Parameter 


Model-A 


Model-B 


Model-C 




0.146 


0.146 


0.146 


is 


7.1 


100.0 


0.1 ~7~ 50 


T 


0.74 


0.05 


0.10 


tgu, 


5.12 


0.39 


3.58 


Z-max 


0.0964 


0.0713 


0.0878 


(Z) 


0.0365 


0.0279 


0.0294 


G(t) 


0.004 


0.002 


0.003 


S(t) 


0.845 


0.942 


0.994 


N en h 


45.7% 


85.2% 


74.8% 



Fig. 4. Model-C: the star formation rate as a function of the 
time for the central region of the galaxy model with 3Ml,i2- 
At the age of 3 Gyr the efficiency of the SFR is let increase 
from v — 0.1 up to v = 50 over a time scale of 10 8 yr. The 
parameters of this model are given in Table 2 




e(Gyr) 



Fig. 5. Model-C: Panel (a) shows maximum and mean metal- 
licity dotted and solid line, respectively. Panel (b) shows the 
fractionary density of gas G(t) and living stars S(t) as a func- 
tion of time, dotted and solid line, respectively 



dances by Anders & Grevesse (1989), Grevesse (1991) and 
Grevesse & Noels (1993), we find the general relation 



'Fe' 


= log 






" a ' 




' a " 








-0.05 


H 






IFe- 


IFei 



(13) 



where the term [a/Fe] stand for all a-elements lumped 
together. 

The recent empirical calibration by Borges et al. (1995) 
for the Mg 2 index includes the effect of different [Mg/Fe] 
ratios 



lnM 52 = -9.037+ 5.795 



0.389 



5040 



T, 



eff 



+ 0.398 log g+ 



'Fe' 


- 0.16 


'Fe' 


2 


\Mg~\ 




+ 0.981 


H 


H 


. Fe _ 



(14) 



which holds for effective temperatures and gravities in the 
ranges 3800 < T eff < 6500 K and 0.7 < \ogg < 4.5. 

To our knowledge, no corresponding calibration for 
the (Fe) index is yet available, so that one is forced to 
adopt that with no dependence on [Mg/Fe). Neverthe- 
less, a zero-order evaluation of the effect of [Mg/Fe] on 
the (Fe) index is possible via the different relation be- 
tween Z and [Fe/H] of the [a/Fe] + 1 case. 

The above relations for [Fe/H], Mg 2 and implicitly 
(Fe) are used to generate new SSPs and galactic models 
in which not only the chemical abundances are enhanced 
with respect to the solar value but also the effect of this 
on the line strength indices is taken into account in a self- 
consistent manner. 
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4-2. Three different chemical structures 

We present here three galactic models that in virtue of 
their particular history of star formation, have differ- 
ent chemical structures and degree of enhancement in 
a-elements. The discussion is limited to the central region 
of the 3Mt 12 galaxy. 

Model- A: late galactic wind. This case has late galactic 
wind (~ 5.12 Gyr), which means that Type la supernovae 
dominate the enrichment in Fe of the gas, and the ratio 
[a/Fe] is solar or below solar for most of the time. This 
model is actually the central region of the 3 Ml. 12 galaxy 
presented above. The percentage (N en h) of a-enhanced 
stars that are still alive at the age of 15 Gyr amounts to 
45.7%. 

Model-B: early galactic wind. In order to enhance the 
relative abundance of elements from Type II supernovae 
we arbitrarily shortened the duration of the star form- 
ing period. To this aim, in the central region of the same 
galaxy, the efficiency of star formation has been increased 
v = 100) and the infall time scale decreases (r=0.05 Gyr) 
so that the galactic wind occurs much earlier than in the 
previous case (at 0.39 Gyr). The material (gas and stars) 
of Model-B is therefore strongly enhanced in a-elements. 
The percentage (N en h) of a-enhanced stars that are still 
alive at the age of 15 Gyr amounts to 87.2%. 

Model- C: recent burst of star formation. A third possi- 
bility is considered, in which a burst of star formation can 
occur within a galaxy that already underwent significant 
stellar activity and metal enrichment during its previous 
history. This model (always limited to the central region of 
the galaxy) has nearly constant star formation rate from 
the beginning, but at the age (arbitrarily chosen) of 3 
Gyr it is supposed to undergo a sudden increase in the 
star formation rate. To this aim, the specific efficiency of 
star formation v is let increase from v = 0.1 to v = 50 
over a time scale of 10 8 yr. The initial nearly constant 
stellar activity is secured by adopting a long time scale of 
gas accretion in the infall scheme (r = 10 Gyr). The rate 
of star formation (in units of MQyr -1 ) as a function of 
time is shown in Fig. ^. Soon after the intense period of 
star formation, the galactic wind occurs thus halting star 
formation and chemical enrichment. The basic chemical 
properties of the model as a function of the age are shown 
in Fig. g. This displays the maximum (dotted line) and 
mean metallicity (solid line), and the fractionary mass of 
gas G(t) (dotted line) and stars S(t) (solid line) both nor- 
malized to AMl^s- The percentage (N en h) of a-enhanced 
stars that are still alive at the age of 15 Gyr amounts to 
74.5%. 

The basic data for the three models in question are 
summarized in Table 2, whereas the evolution of their 
chemical abundances and present-day partition function 
N(Z) are shown in Figs. |[ where the left panel is for 



Model- A, the middle panel is for Model-B, and the right 
panel is for Model-C. These figures are the analog of Fig. |[ 

4.8. Mg 2 and (Fe) in SSPs with [a/Fe] ^ 

To secure full consistency between chemical abundances 
and line strength indices for the galactic models and make 
use of the SSP technique (cf. Bressan et al. 1996) one 
should adopt SSPs with the same pattern of abundances 
indicated by the chemical models. 

To this aim, the distributions of chemical abundances 
(and their ratios) as a function of the total metallicity 
and/or time provided by the model galaxies are used to 
calculate the Mg 2 and (Fe) indices of SSPs with different 
Z, [Fe/H], [O/Fe], and [Mg/Fe]. 
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Fig. 7. Panels (a), (c) and e) show the Mg2 index evolution for 
SSPs with different metallicity (Z=0.0004, Z=0.004, Z=0.008, 
Z=0.02, Z=0.05 and Z=0.1 solid dotted dashed long-dashed, 
dot-dashed and dot-long-dashed lines, respectively) and the as- 
sumption of enhancement in a-elements. Panels (b), (d) and 
(f) show the same but without enhancement of a-elements 

To quantify the abundance of a-elements we prefer to 
use [O/Fe] instead of [Mg/Fe], because the stellar yields 
by Portinari et al. (1997) that are at the base of our chem- 
ical models somewhat overestimate the production of Fe 
by Type II supernovae, and underestimate in turn the ra- 
tio [Mg/Fe] as compared to the observational value. Ac- 
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Fig. 8. Panels (a), (c) and e ) show the (-Fe) index evolution for 
SSPs with different metallicity (Z=0.0004, Z=0.004, Z=0.008, 
Z=0.02, Z=0.05 and Z=0.1 solid dotted dashed long-dashed, 
dot-dashed and dot-long-dashed lines, respectively) and the as- 
sumption of enhancement in ct-elements. Panels (b), (d) and 
(f) show the same but without enhancement of a-elements 



cording to Portinari et al. (1997) the theoretical [Mg/Fe] 
is about 0.2-^0.3 dex lower than indicated by the observa- 
tional data. Applying this correction, the ratio [Mg/Fe] 
gets close to the ratio [O/Fe] (see Fig. which some- 
how justifies our use of [O/Fe] instead of [Mg/Fe] in the 
procedure below. This marginal drawback of the chemi- 
cal model does not however affect the conclusion of this 
analysis. 

For any value of Z we read from Fig. ^| the ratios 
[C/FeL [O/Fe], and [Mg/Fe], derive [Fe/H] from equa- 
tion (§J), and insert [Mg/Fe], i.e. = [O/Fe], and [Fe/H] 
into equation (jl4|). It goes without saying that the val- 
ues of [Fe/H] derived from Figs. || and equation ( |l3|) are 
mutually consistent by definition. 

Table 3 shows the values of [Fe / H] and [O / Fe] as- 
signed to the SSPs according to the chemical structure 
of Model-A, Model-B, and Model-C and, for purposes of 
comparison, to reference SSPs with no enhancement at all. 

The temporal evolution of the Mg 2 and (Fe) indices 
for the SSPs listed in Table 3 is shown in the various panels 
of Figs. and ||, respectively. 



It is soon evident that with [a/Fe] = (right panels of 
Fig. Q), Mg 2 monotonically increases with the metallicity, 
but for the extreme SSP with Z=0.1 for which the trend 
is reversed at ages older than 5 Gyr. When [a/Fe] ^ 0, 
the trend is more complicated as it depends on the degree 
of enhancement. For cases A and B, the strongest Mg2 
happens to occur for Z = 0.02, whereas for case C it 
occurs for Z = 0.05 (left panels of Fig. [?]). 

As far as (Fe) is concerned, with [ct/Fe] = the index 
gets stronger at increasing metallicity but for the extreme 
case of Z = 0.1, in which (Fe) gets lower than or com- 
parable to the values for the case with Z — 0.05 at ages 
older than about 5 Gyr (right panels of Fig. |J). In pres- 
ence of enhancement in a-elements, there is a significant 
dependence of (Fe) on this parameter even at our zero- 
order evaluation (see the left panels of Fig. ||). Although 
the exact evaluation of the effect of [Mg/Fe] on (Fe) is 
hampered by the lack of the proper calibration, still the 
above experiments clarify that it cannot be neglected. 

4--4- The Mg2 and (Fe) indices of galaxies 

What the results for the Mg2 and (Fe) indices would be 
when applying these SSPs to model galaxies ? The situ- 
ation is displayed in the various panels of Fig. || for Mg 2 
and Fig. [H^ for (Fe) . The combined analysis of the chem- 
ical structures, partition functions N(Z), and temporal 
variations of the Mg 2 and (Fe) indices of the model galax- 
ies (Figs. ^ and [H]) allow us to make the following re- 
marks: 

(i) [Mg/Fe] ^ 0: Mg 2 in Model-A (late wind, no chem- 
ical enhancement) is always weaker than in Modcl-B 
(early wind, significant chemical enhancement). How- 
ever, the difference is large for ages younger than about 
5 Gyr, and gets very small up to vanishing for older 
ages. Mg 2 of Model-C is always weaker than Model-A 
and Model-B which means that the higher (or compa- 
rable) enhancement in a-elements of Model-C with re- 
spect to the previous ones does not produce a stronger 
Mg 2 index. The age dependence of (Fe) is more in- 
trigued. At young ages (< 5 Gyr) the intensity of (Fe) 
gets stronger passing from Model-C to Model-A and 
finally Model-B. At ages older than 5 Gyr, Model-A 
has the strongest (Fe) whereas Model-B and Model-C 
are weaker and of the same intensity. 

(ii) [Mg/Fe] — 0: Mg 2 in Model-A (late wind, no chem- 
ical enhancement) is first weaker than in Model-B 
(early wind, significant chemical enhancement) up to 
ages of about 5 Gyr, and then becomes significantly 
stronger at older ages. Finally, the Mg 2 index of 
Model-C (burst of star formation, and strong enhance- 
ment) is weaker or about equal to that of Model-A and 
Model-B. The index (Fe) closely follows the trend of 
Mg 2 all over the age range. 

It follows that both Mg 2 and (Fe) don not simply cor- 
relate with age, [Fe/H] and [Mg/Fe]. The striking result 
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Table 3. [O/Fe] and [Fe/H] ratios for SSPs with enhancement of a-elements according to Model-A, Model-B and Model-C. 
The same ratios for the reference SSPs with no enhancement are also shown 



z 


Model-A 


Model-B 


Modol-C 


Reference SSPs 




[O/Fe] 


[Fe/H] 


[O/Fe] 


[Fe/H] 


[O/Fe] 


[Fe/H] 


[O/Fe] 


[Fe/H] 


0.0004 


+0.8 


-2.38 


+3.28 


-4.87 


+0.51 


-2.13 


0.0 


-1.71 


0.004 


+0.6 


-1.22 


+ 1.72 


-2.23 


+0.14 


-0.82 


0.0 


-0.71 


0.008 


+0.5 


-0.80 


+ 1.22 


-1.44 


+0.03 


-0.42 


0.0 


-0.39 


0.02 


+0.4 


-0.30 


+0.72 


-0.57 


+0.30 


-0.22 


0.0 


+0.03 


0.05 


-0.50 


+0.88 


-0.25 


+0.69 


+0.31 


+0.24 


0.0 


+0.50 


0.1 


-0.80 


+ 1.55 


-0.59 


+ 1.40 


-0.87 


+ 1.60 


0.0 


+0.95 



is that Mg 2 and to some extent (Fe) as well of mod- 
els with supposedly the highest degree of enhancement in 
a-elements happen to be weaker than in those with no 
enhancement. 

What causes the odd behaviour of Mg 2 and (Fe) as a 
function of the age and underlying chemical structure of 
the model galaxy ? 

To answer this question, we have artificially removed 
from the partition function N(Z) all the stars in certain 
metallicity bins and re-calculated the line strength indices 
for the three models. 

Panels (c) and (d) of Figs. |j] and [h] show the results 
when all stars with metallicity higher than Z=0.05 are 
removed. This is motivated by the trend as function of 
the metallicity shown by the SSPs we have already pointed 
out. The situation remains practically unchanged. 

Likewise, panels (e) and (f) of Figs. ^|and show the 
same but when all stars with metallicity lower than Z = 
0.008 are removed. Now the results change significantly. 
In the case of [Mg/Fe] ^ 0, the Mg 2 index of Model-B 
is always much stronger than that of Model-A. In such 
a case there is correspondence between the strength of 
the index and the amount of enhancement in et-elements. 
In contrast, (Fe) of Model-A is first slightly weaker and 
then stronger than in Model-B, whereas that of Model-C 
is first weaker and then almost equal to that of Model- 
B. cases. In the case of [Mg/Fe] — 0, at ages older than 
about 5 Gyr the Mg 2 and (Fe) of the three models are 
almost equal each other, with a marginal increase passing 
from Model-C to Model-B and Model-A. At younger ages, 
while Model-A and Model-C have nearly the same indices, 
Model-B has always the strongest values. 

The above analysis clarifies in a quantitative fashion a 
number of important clues: 

— In galaxies, both Mg 2 and (Fe) depend on the age, 
N(Z), [Fe/H], and [Mg/Fe] in a somewhat unpre- 
dictable fashion. 

— Because of this, inferring from the above indices the 
abundance ratio [Mg/Fe] and the enhancement of a- 



elements is a difficult task with somewhat ambiguous 
results. 

— Different slopes of the gradients in Mg 2 and (Fe) do 
not automatically imply gradients in chemical abun- 
dances or enhancement ratios. 



Mg, ([Mg/Fe] - 0) (Fe) ([Mg/Fe] » 0) 




0.0002 0.002 0.02 0.0002 0.002 0.02 

z z 



Fig. 11. Curves of constant Mg2 (left) and (Fe) (right) at 
varying age (in Gyr) and metallicity (Z) for SSP with no en- 
hancement in a-elements. The calibrations in use are from 
Worthey et al. (1994). The intensity of the indices is anno- 
tated along the curves 

5. Indices in the age-metallicity plane of SSPs 

We like to conclude this study presenting the loci of con- 
stant Mg 2 and (Fe) in the age-metallicity plane of SSPs. 
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Mg, ([Mg/Fe] * 0) (Fe) ([Mg/Fe] * 0) 




Fig. 12. The same as in Fig. 11 but for SSP with enhancement 
if a-elements according to the entries of Table 3. The calibra- 
tions in use is from Borges et al. (1995) for Mg2 and Worthey 
et al. (1994) for {Fe) however corrected for the different re- 
lation between Z and [Fe/H] in presence of enhancement in 
a-elements. The intensity of the indices is annotated along the 
curves 

This is shown in Figs. |ll| and |l2] displaying the above in- 
dices with and without enhancement of a-elements using 
the calibrations discussed in the previous sections. 

This plane can be used to quickly check how gradients 
in ages and/or metallicities across galaxies (inferred from 
other independent analysis) would reflect onto gradients 
in Aft/2 and (-Fe) and vice-versa. 

An interesting feature to note, is the marked dip in 
the Mg 2 index for values greater than a certain limit that 
occurs at a certain value of the metallicity. The thresh- 
old value for Mg 2 is about 0.26 and the metallicity is 
about 0.03 in the case of no enhancement, whereas they 
are lowered to 0.15 and 0.008, respectively, in presence of 
enhancement. 

Starting from the observational result that in gen- 
eral elliptical galaxies show stronger Mg 2 and (Fe) to- 
ward the center, and representing the local mix of stellar 
population in a galaxy (center and/or external regions) 
with a mean SSP of suitable age and composition, we see 
that the observational gradient in Mg 2 and (Fe) could 
be compatible with (i) either a nucleus more metal-rich 
and older than the external regions; (ii) or nucleus more 
metal-rich and younger than the external regions. In con- 
trast a nucleus less metal-rich and older than the external 
regions would lead to gradients in Mgi and (Fe) oppo- 



site to what observed. The situation is straightforward 
with the Worthey (1992) calibrations, whereas it is some- 
what intrigued with the Borges et al. (1995) calibration in 
presence of a-enhancement. It worth recalling here that 
Bressan et al. (1996) analyzing the Gonzales (1993) Hp 
and [MgFe] data for elliptical galaxies and their variation 
across these systems suggested that most galaxies ought 
to have nuclei with higher metallicities and longer dura- 
tions of star forming activity than the peripheral regions. 
This suggestion is fully compatible with the gradients in 
Mg2 and (Fe) observed in elliptical galaxies. 

6. Summary and conclusions 

Aim of this study was to ascertain on a quantitative ba- 
sis the effect of age, metallicity, partition function N(Z), 
abundance ratio [Mg/Fe], and calibration in usage on the 
line strength indices Mg 2 and (Fe) from whose gradients 
the problem of the possible enhancement of a-elements to- 
ward the center of elliptical galaxies originates. Although 
by no means we want to exclude such a possibility, the 
attention is called on a number of indirect effects that 
could invalidate the one-to-one correlation between the 
index intensity and the abundance of the corresponding 
element, and indirectly on the one-to-one correlation be- 
tween the relative slopes of the observational gradients and 
the inferred spatial variation of abundance ratios ([ce/Fe] 
in particular). 

The results of this study can be summarized as follows: 

1. The intensity of Mg 2 does not simply correlate with 
the abundance of Mg and the ratio [Mg/Fe] in par- 
ticular. 

2. The intensity of Mg 2 does not simply correlate with 
the age or the metallicity. 

3. The intensity of Mg 2 much depends on the partition 
function N(Z). 

4. Likewise for the (Fe) index. 

5. Inferring the abundance of Mg or the enhancement 
ratio [Mg / Fe] is a cumbersome affair whose solution is 
not always possible because hints on N(Z) are needed. 

6. The observational gradients in Mg 2 and (Fe) do not 
automatically imply gradients in the abundances of 
Mg and Fe or enhancement ratios. Inferring from the 
observational Mg 2 and (Fe) constraints on the past 
history of star formation (via the different time scales 
of Mg and Fe enrichment) may be risky. 

Although most of these conclusions have already been 
around in literature, their quantitative assessment has 
never been attempted before in a systematic fashion, in 
particular in the complex but realistic situation in which 
many stellar populations with different ages, metallicities, 
and abundance ratios are present. 
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Fig. 2. Gradients in colors and line strength indices for the galaxy NGC 6407 (Carollo & Danziger 1994a). Left Panel: comparison 
with the theoretical gradients in (B-R) for models of different mass and same age (15 Gyr). Right Panel: comparison with the 
theoretical gradients in line strength indices for the 3 x Ml, 12 model which has nearly the same M/Lb ratio as NGC 6407. 
Finally, the observational data along major and minor axis are indicated by full and empty circles, respectively 
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Fig. 3. Panels (a) and (b): the number of living stars and abundance ratio distribution per metallicity bin in the central core 
and the I s * shells, respectively, of the galaxy with mass 3Ml,i2- The abundance ratios are in the standard notation. The solid 
line is distribution of living stars in units of 1O 11 M0. The dotted, dashed, long-dashed, and dot-dashed lines give the distribution 
per metallicity bin of [O/Fe], [Fe/H], [Mg/Fe], and [C /Fe], respectively. The top scale gives the birth-time t = To — Tssp in 
Gyr of a SSP of age Tssp in a galaxy of age Tq 
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Fig. 6. The partition function N(Z) and abundance ratios distribution per metallicity bin, for the Model- A (left panel), Model-B 
(middle panel), and Model-C (right panel). The solid line is N(Z) in units of 10 11 at the age of 15 Gyr. The dotted, dashed, 
long-dashed, and dot-dashed lines give the distribution per metallicity bin for [O/Fe], [Fe/H], [Mg/Fe] and [C/Fe], respectively. 
The abundance ratios are in the standard notation. The top scale gives the birth-time t — To — Tssp of a SSP with age Tssp 
in a galaxy with age Tg 




Fig. 9. Evolution of the Mg2 index as a function of time. Panels (a), (c) and (e) show the evolution of the Mg? index 
calculated including the effect of the chemical abundances, while Panels (b), (d) and (f) show the same but without the effect 
of a-enhancement. The solid line corresponds to Model-A the dotted line to Model-B, and dashed line to Model-C 
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Fig. 10. Evolution of the (Fe) index as a function of time. Panels (a), (c) and (e) show the evolution of the Mg2 index 
calculated including the effect of the chemical abundances, while Panels (b), (d) and (f) show the same but without the effect 
of a-enhancement. The solid line corresponds to Model-A the dotted line to Model-B and dashed line to Model-C 



